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Relaxation of writhe and twist of a bi-helical magnetic field 
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Abstract. In the past few years suggestions have emerged that the solar magnetic field might have a bi-helical contribution 
with oppositely polarized magnetic fields at large and small scales, and that the shedding of such fields may be crucial for the 
operation of the dynamo. It is shown that, if a bi-helical field is shed into the solar wind, positive and negative contributions 
of the magnetic helicity spectrum tend to mix and decay. Even in the absence of turbulence, mixing and decay can occur on a 
time scale faster than the resistive one provided the two signs of magnetic helicity originate from a single tube. In the presence 
of turbulence, positively and negatively polarized contributions mix rapidly in such a way that the ratio of magnetic helicity 
to magnetic energy is largest both at the largest scale and in the dissipation range. In absolute units the small scale excess of 
helical fields is however negligible. 
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1. Introduction 

Hydromagnetic turbulence in the solar wind is long known 
to have a magnetic power spectrum that is compatible with 
(although with substantial error bars, see Matthaeus et 
al. 1982). The magnitude of the magnetic helicity is small, 
but there is some evidence that it is negative above the helio- 
spheric current sheet and positive below it; see Smith & Bieber 

(1993) and Bieber & Rust (1995). Spectra of magnetic energy 
and helicity have later also been calculated by Goldstein et al. 

(1994) , Leamon et al. (1998), and Smith (1999) under the as- 
sumption of isotropy. They found net magnetic helicity only in 
the dissipation range. They also reported efficient cancellation 
of magnetic helicity in the inertial range. 

The presence of magnetic helicity at small scales is puz- 
zling, because magnetic helicity is known to cascade to large 
scales (Frisch et al. 1975, Pouquet et al. 1976). This inverse 
cascade has also been seen in turbulence simulations both with 
helical forcing (Pouquet & Patterson 1978, Balsara & Pouquet 
1999, Brandenburg 2001) and without forcing but helical ini- 
tial fields that are either tube-like (Horiuchi & Sato 1985, 1986, 
Zhu et al. 1995) or turbulent (Christensson et al. 2001). It is 
possible, however, that the spectral transfer of magnetic energy 
from small to large scales is a process that does not happen 
on a dynamical time scale but on a much slower resistive time 
scale (Brandenburg & Sarson 2002). In that case the small scale 
magnetic helicity will not have had time to cascade to larger 
scales in the interval between the time of ejection of the ob- 
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served magnetic field patch from the solar surface and the time 
of in-situ observation. 

The purpose of the present paper is to study the type of in- 
teractions that can occur on short dynamical time scales. This 
is important in connection with the solar magnetic field. In ac- 
tive regions the field is found to be markedly helical (Seehafer 
1990, Pevtsov et al. 1995, Rust & Kumar 1996, Bao et al. 
1999). On the other hand, the observed absence of net mag- 
netic helicity in the inertial range of the solar wind seems to be 
in conflict with the observed markedly helical field on the solar 
surface. We believe that this may be related to recent theoreti- 
cal suggestions that the magnetic helicity ejected from the sun 
is bi-helical and has two components with opposite magnetic 
helicity: one from intermediate scale fields and one from the 
global field of the sun (Blackman & Brandenburg 2003). 

The bi-helical nature of the solar magnetic field is high- 
lighted by the fact that in active regions the magnetic he- 
licity is negative (and positive on the southern hemisphere) 
while, on the other hand, bipolar regions are tilted according 
to Joy's law (Hale et al. 1919), i.e. in the clockwise direction 
in the north (and anti-clockwise in the south) suggesting that 
the swirl of flux tubes is right handed, giving rise to positive 
magnetic helicity (and negative magnetic helicity in the south). 
This corresponds to what is also known as writhe helicity (e.g. 
Longcope & Klapper 1997, Demoulin et al. 2002). This idea 
of a bi-helical field is further supported by studies of sigmoids. 
Figure 2a of Gibson et al. (2002) shows a TRACE image of 
an N-shaped sigmoid (right-handed writhe) with left-handed 
twisted filaments of the active region NOAA AR 8668, typical 
of the northern hemisphere. 
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The puzzle of a bi-helical field structure can be resolved 
when one realizes that an overall tilting of flux tubes causes 
simultaneously an equal amount of internal twist in the tube 
(Blackman & Brandenburg 2003). This is best seen when the 
magnetic field is pictured as a ribbon. This picture supports the 
notion that the magnetic field generated in the sun by an a- 
effect must have opposite signs of magnetic helicity at large 
and small scales (Seehafer 1996, Ji 1999, Brandenburg 2001, 
Field & Blackman 2002, Blackman & Brandenburg 2002). 

Once a magnetic field structure has emerged at the solar 
surface, the tilt angle is known to relax gradually on a time 
scale of a few days (Howard 1996, Longcope & Choudhuri 
2002). This means that also the internal twist must decreases. In 
order to understand quantitatively the fate of a bi-helical mag- 
netic field after it has been ejected into the solar wind we study, 
using simulations of hydromagnetic turbulence, what happens 
to a magnetic field that is composed of large and small scale 
fields of opposite helicity and different amplitudes. The expan- 
sion of the solar wind is not explicitly taken into account, but it 
is known that under similar circumstances of cosmological ex- 
pansion helical magnetic fields can still show the inverse cas- 
cade effect (Brandenburg et al. 1996). 

2. The model 

We consider a compressible isothermal gas with constant sound 
speed c s , constant kinematic viscosity v, constant magnetic dif- 
fusivity rj, and constant magnetic permeability /.to- The govern- 
ing equations for density p, velocity u, and magnetic vector 
potential A, are given by 



Dlnp 
Dt 



= -V • M, 



Du 9 J x B 

— = -c s 2 Vlnp+ 

Dt p 



-^visc ~t~ f 



dA 
~dt 



11 x B — rjp J, 



(1) 

(2) 
(3) 



where D/Dt = d/dt + u ■ V is the advective derivative, B — 
V x A is the magnetic field, and J = V x B/ po is the current 
density. The viscous force is 



-FVisc = v (V 2 m+ ±VV • u + 2^S • Vlnp) 



(4) 



i<5y V • u is the traceless rate 



where Sy = \{uij + Uj.i 
of strain tensor. In cases with forcing, / ^ is a nonhelical 
forcing function, selected randomly at each time step from a 
set of nonhelical transversal wave vectors 



/ fe = (fc x e) A/* 2 - (*: • e) 2 



(5) 



where e is an arbitrary unit vector needed in order to generate 
a vector k x e that is perpendicular to k. 

We use nondimensional quantities by measuring u in units 
of c s , x in units of 1/ki, where fci is the smallest wave number 
in the box (side length L — 2tt), density in units of the ini- 
tial value po, and B is measured in units of ^popo c s . This is 
equivalent to putting c s = ki = po = A*o = !• 



In a periodic domain the total helicity, 



H = A ■ B dV, 



(6) 



is gauge invariant and conserved in the limit of zero magnetic 
diffusivity. H is a topological measure of the mutual linkage of 
the magnetic flux lines and thus the complexity of the magnetic 
field. 

Since we are interested in the distribution of helicity be- 
tween different scales it is convenient to decompose the total 
magnetic helicity into spectral modes, Hk, normalized such 
that 



H = H k dfc. 



(7) 



Hk is the contribution to H from the wave number interval (e.g. 
Brandenburg et al. 2002) 



k-5k/2< |fc| < k + Sk/2 (fc-shell). 



3. Results 



(8) 



3.1. No forcing 

We first consider the evolution of a simple bi-helical initial field 
consisting of a superposition of two Beltrami waves with wave 
numbers k\ = (0, 0, 1) and k§ = (5, 0, 0) and helicities Hi = 
10~ 6 and H§ = — 10 -5 . The vector potential is then 








cos k$x 
' sinkc,x 



(9) 



The net magnetic helicity is negative, but since the small scales 
decay faster there will be a time when the net magnetic helicity 
turns positive. 

Figure[2shows magnetic and kinetic power spectra, 



E k = h J \ u k\ 2 ^ M * = \ J \ B k\ 



■dk, 



(10) 



fc-shell 



fc-shell 



at three different times. We find (as expected) resistive decay 
of the power on the two wave numbers. Thus, there is no en- 
hanced decay. This is readily explained by the absence of any 
spectral overlap between the two components. This would not 
change even if the wave vectors of the two Beltrami waves were 
pointing in the same direction. 

3.2. Nonhelical forcing 

We now drive turbulence by adding a random forcing term that 
acts on large scales with wave numbers between 1 and 2. We 
find that the magnetic power in the two modes at fc = 1 and 
5 spreads rapidly among other wave numbers; see Fig. In 
addition, there is also dynamo action giving rise to a magnetic 
power spectrum that rises with t; see the last panel of Fig. [2] 
The form of the magnetic energy spectrum in the dynamo case 
is subject of a separate investigation (Haugen et al. 2003, and 
references therein). 
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Fig. 1. Magnetic and kinetic power spectra at three differ- 
ent times. No explicit forcing of the flow. The initial mag- 
netic field consists of a superposition of two Beltrami fields 

at k = (0,0, 1) with H x = 10~ 6 and at k = (5,0,0) with 
H 5 = -10~ 5 . 
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Fig. 2. Magnetic and kinetic power spectra at three different 
times. The flow is forced nonhelically in the wave numbers 
interval 1 < |fc| < 2. The initial magnetic field consists of 
a superposition of two Beltrami fields at k = (0,0, 1) with 
Hi = 10~ 6 and at k = (5, 0, 0) with H 5 = -10~ 5 . 

It is convenient to divide the magnetic power spectrum 
into contributions from positively and negatively polarized 
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Fig. 3. Magnetic power spectra of positive and negative polar- 
ized contributions to the magnetic field. The sign of the mag- 
netic helicity is indicated by black dots (positive) and open 
squares (negative). Otherwise like in Fig. [2] 

waves such that M k = + (Waleffe 1993, see also 
Brandenburg et al. 2002). The magnetic helicity spectrum can 
then be written as 



H k = (2/k)(M+-M^). 



(11) 



The result is shown in Fig. [3] where we see that at early times 
(t < 10), peaks at k = 1 and MjT peaks at k = 5. The 
magnetic helicity spectrum, normalized by k /2, is positive for 
k < 3 and negative for k around 5. At later times, sa MZ 
and so the normalized magnetic helicity, ^kHk, is small, sug- 
gesting that contributions from positively and negatively polar- 
ized waves have mixed almost completely. 

In Fig.|4]we plot ^kHk/Mh, i.e. the magnetic helicity spec- 
trum relative to the magnetic energy spectrum. Note the sys- 
tematically negative values of the normalized magnetic helicity 
in the dissipative subrange, just like in the plots of Leamon et 
al. (1998). This enhancement of (negative) magnetic helicity is 
probably the result of a direct cascade that transfers the nega- 
tive magnetic helicity from intermediate to still smaller scales. 
On the other hand, in absolute units the magnetic helicity re- 
mains small at small scales and therefore the enhancement rel- 
ative to the magnetic energy spectrum could be regarded as a 
typical feature of such normalized plots and is not an indication 
of a real excess of small scale magnetic helicity. 

To summarize, the decay of magnetic helicity by mixing 
occurs on a resistive time scale if there is no externally driven 
turbulence. Decay on a faster (turbulent) time scale is possible 
in the presence of externally driven nonhelical turbulence. This 
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Fig. 4. Magnetic helicity spectrum normalized by magnetic en- 
ergy for the same run as shown in Fig. [2] at t = 100. Note the 
relative magnetic helicity excess both at largest and smallest 
scales. 



is also in qualitative agreement with simulations by Maron & 
Blackman (2002) who found that if the helicity of the forcing is 
below a certain threshold the large scale field decays rapidly. In 
the remainder of this paper we consider the decay of bi-helical 
fields without any driving. 



3.3. Decay of a random bi-helical field 

We first consider a random initial condition obtained by 
stopping a helically driven turbulence run (as studied in 
Brandenburg 2001). In Fig.[5]we show the evolution of Mt and 
M^T for those values of k where the spectra are maximum for 
the initial condition (k = 1 for and k = 5 for M^). Note 
that both Mfc and MjT decay at rates that are clearly slower 
than the resistive rate (indicated by a straight line). This is be- 
cause in three dimensional turbulence the decay of a magnetic 
field is easily offset by dynamo action, even if the velocity is 
decaying (Dobler et al. 2003). 

In addition to the field decaying only slowly we also find 
that the bi-helical character of the field prevails. Thus, dynamo 
action alone does not produce mixing of the positively and neg- 
atively polarized components. In the following section we sug- 
gest that this is primarily a shortcoming of our initial condition 
lacking positively and negatively polarized components within 
a single physically connected flux structure. 



3.4. Decay of a bi-helical flux tube 

In order to make the connection with the magnetic field struc- 
ture in real space it is instructive to consider explicitly a twisted 
flux tube. We do this here by constructing the Cauchy solution 
of an initially straight flux tube in a simple steady flow field. 
The Cauchy solution (e.g. Moffatt 1978) is 
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Fig. 5. Decay of (t) f or k = 1 (dotted lines) and k = 5 
(solid lines). The lower curves, denoted by (i), are for v = 
2 x i) = 10~ 3 while the upper ones, denoted by (ii), are for 
v = r\ = 10~ 3 . The run that was forced helically until t = 5, 
and was then let to decay. Note that the field never decays faster 
than on the resistive rate (here only indicated for k — 5). 




Fig. 6. Magnetic flux tube constructed from the Cauchy solu- 
tion, t = 2, ip = 0.2. 

where r is the time coordinate for constructing the solu- 
tion, Gij = dxi /dxoj is the Lagrangian displacement matrix, 
Bq(xq) = B(x, 0) is the initial condition, and 



x(x q ,t) =x Q + [ U {x{t'))At' 
Jo 



(13) 



is the position of an advected test particle whose original posi- 
tion was at xq. 

The flow must have the property of lifting the tube only in 
one section and tilting it sideways, hence causing writhe helic- 
ity in the tube. A simple steady flow with such properties is 





U(x) — cpz sinx 
1 + cos x 



(14) 



B 1 (x,t) 



Gjj(x ,T) 

detG 



B j(x ), 



(12) 



where <p is a parameter that controls the amount of twist. 
Periodic boundary conditions are assumed in the range — 7r < 
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Fig. 7. Initial spectra of magnetic helicity, Hf., and of magnetic 
energy of positively and negatively polarized components, 
and Mfr ', respectively. 




Fig. 8. Magnetic flux tube after letting it relax until t = 33. 

(x, y) < it and < z < 2tt. With this velocity field, Eq. fljl 
reduces to 

x(x ,t) =x + U(x )t. (15) 
This yields 

/ 1 \ 

Gy(xo)T")= zoiprcosxo 1 iprsinxo . (16) 
\ — Tsinxo 1 J 

Note that dct G = 1 (incompressibility) and Gij(xa, 0) = Sy, 
so Bo(xo) = B(x,0). An example of the resulting magnetic 
field structure is shown in Fig.|6]for ip = 0.2 and r = 2. The 
initial field was a straight tube with 

B {x) = xB Q exp{-[y 2 + (z ~ h) 2 ]/d 2 }, (17) 

where ft, = 1 is the initial height of the tube, d = 0.5 is its 
radius, and Bo = 0.1 is the initial field strength. The magnetic 
helicity spectrum shows distinctively bi-helical behavior; see 

Fig.m 

We use this Cauchy solution as initial condition for our sim- 
ulation, i.e. t = corresponds now to r = 2. The result at time 
t = 33 is shown in Fig. [8] Unlike the previous case (Figs[2and 
0, the magnetic energies of the positively and negatively po- 
larized components decay at a rate faster than the resistive rate. 
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Fig. 9. Decay of Mf (t) and M/(t) for v = r} = 10" 3 
(solid lines). The dotted lines give the result for v — 10 -3 and 
T) = 10~ 4 . The slopes (dash-dotted lines) are given in units 
of 7/ t = Mrms/fceff, where M rms is the peak rms velocity and 
fc e ff = 4 is the estimated effective wave number of the energy 
carrying scale during the time where the velocity is maximum. 
The upturn in the M^(t) curves after t = 40 is temporary and 
caused by field cascading backward from larger wave numbers. 

We characterize the decay of M± (k±,t) at fc_ = 1 and k + = 4 
in terms of effective diffusivity coefficients, r]±, by writing the 
decay law as 

M ± (k±,t) - exp(-2r?±fe|t). (18) 

The result is shown in Fig. [9] for two different values of ij. 
The decay rates are either completely independent of the mi- 
croscopic value of r\ (for k = 1), or depend on it only weakly 
(for k = 4). The decay happens therefore on something like a 
'turbulent' resistive time scale. For orientation we use the max- 
imum rms velocity u rms sa 0.006 and the estimated effective 
wave number of the energy carrying scale during velocity max- 
imum, fc ff = 4, to calculate a 'turbulent' magnetic diffusivity, 
rjt = Uj-ms/fceff. It turns out that the effective diffusivity co- 
efficients, r)±, are of the order of rj t and independent of the 
magnetic Reynolds number. A relevant measure of the mag- 
netic Reynolds is rit/rj. This ratio is either 1.5 or 15 in the two 
cases shown (Fig. [9] solid and dashed lines, respectively), so 
the magnetic Reynolds number has changed by a factor of 10 
while the actual decay rates (characterized by the values of r]+ 
and ?7_) are almost unchanged (Fig. |9j. 

We conclude from this that a single flux tube in a bi-helical 
state with writhe and twist helicity of opposite sign relaxes to 
a nonhelical state on a dynamical time scale. The final state is 
not, however, a simple straight tube, but one with a more com- 
plicated internal structure. On the other hand, bi-helical fields 
generated randomly decay more slowly and do not tend to mix. 
This is explained by dynamo action of the decaying turbulent 
flow field and by inverse transfer of magnetic helicity (Pouquet 
et al. 1976, Christensson et al. 2001). 

4. Discussion 

Helical driving in the solar convection zone tends to produce 
bi-helical fields. In the northern hemisphere the magnetic he- 
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licity is probably positive at large scales and negative at smaller 
scales. In the absence of helical driving, for example in the so- 
lar wind, an initially bi-helical field can relax in such a way that 
contributions with positive and negative magnetic helicity mix. 
The decay of magnetic helicity in the different modes will in 
general happen on a turbulent time scale, until a certain level 
is reached where a nonhelical magnetic field can be sustained 
by nonhelical dynamo action. Even in the absence of any driv- 
ing magnetic helicity of opposite signs and at different scales 
can mix, provided this involves a singly connected magnetic 
field structure (as shown in Sect. 13 .4> . This is probably related 
to the possibility of rapid propagation of twist along the tube 
as a torsional Alfven wave (Longcope & Klapper 1997). If the 
different signs of helicity do not reside on a single tube, the 
mixing of oppositely oriented writhe and twist can be as slow 
as a resistive time scale, as is the case where two orthogonal 
Beltrami waves with opposite sign and different wavelengths 
are superimposed. 

In order to see the distribution of the sign of the magnetic 
helicity over different scales one often plots the magnetic he- 
licity spectrum normalized by the magnetic energy spectrum. 
For the solar wind above the heliospheric current sheet this ra- 
tio fluctuates around zero and has a negative excess within the 
dissipative subrange (Leamon et al. 1998). Similar behavior is 
also found in the present case and seems to be indeed a conse- 
quence of magnetic diffusion helping mix oppositely oriented 
waves efficiently. 

Although there is qualitative support of the idea that there 
is large scale magnetic helicity with opposite sign (Gibson et 
al. 2002), there is to our knowledge no quantitative measure- 
ment. All explicit measurements point to a negative sign in the 
northern hemisphere. If the solar field is indeed bi-helical, there 
may be as yet undetected positive magnetic helicity at larger or 
perhaps even at smaller scales. The fact that in the dissipative 
range of solar wind turbulence the magnetic helicity is negative 
supports the former suggestion that the negative helicity seen at 
the solar surface is related to the small scale magnetic helicity, 
and that there is positive magnetic helicity on scales larger than 
the scale of active regions. 

Ideally, one should try to measure this large scale magnetic 
helicity from the longitudinally averaged mean field. From syn- 
optic maps one can obtain the longitudinally averaged radial 
field at the surface. This allows one to calculate the toroidal 
component of the magnetic vector potential, A$. Using spher- 
ical polar coordinates for the axisymmetric mean field, the 
gauge-invariant magnetic helicity of Berger & Field (1984) is 
simply H = 2 J A^B^ dV (see Brandenburg et al. 2002). In 
order to calculate H one still needs B^. Unfortunately, this can- 
not be observed directly. Only its sign can be inferred from 
the orientation of bipolar regions. Preliminary investigations 
(Brandenburg et al. 2003) suggest that there are approximately 
equally long intervals where A^B^ is positive and negative 
during one cycle. The result is therefore inconclusive. 

The best support for oppositely oriented magnetic helic- 
ity at large scales comes from the fact that the footpoints of 
emerging magnetic flux tubes are tilted clockwise in the north- 
ern hemisphere and counterclockwise in the southern hemi- 
sphere (Joy's law). This corresponds to right-handed tubes in 



the northern hemisphere and left-handed tubes in the southern 
hemisphere, i.e. to positive magnetic helicity in the north and 
negative magnetic helicity in the south. This is also in agree- 
ment with studies of sigmoids suggesting that on the northern 
hemisphere N-shaped sigmoids with right-handed writhe have 
left-handed twisted filaments (Gibson et al. 2002). 
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